) 

) 



I 



i 

! 



I 

i 

I 

I 



j 



i' 




% 



A COMPARISON OF FOUR NON-DIVERGENT WIND FIELDS 



by 

Richard Brian Brodehl 

if 

Lieutenant, United States Navy 
B. S. , Naval Academy, 1962 



Submitted in partial fulfillment of the 
requirements for the degree of 

MASTER OF SCIENCE IN METEOROLOGY 

from the 



NAVAL POSTGRADUATE SCHOOL 
September 1967 






J 




( . i 






ABSTRACT 



This paper is concerned with the comparison of four different non- 
divergent wind fields obtained from a single geopotential height field 
over a dense data area. After developing the divergence equation of 
the non-divergent stream function, four different stream functions are 
obtained by modification and/or deletion from the basic equation. 

Isotach patterns for each strea are computed. A comparison 



of the four stream function patterns and their corresponding isotach 
patterns is made. Two of the stream functions are used to obtain the 
non-divergent wind field for the initial wind conditions in a primitive 
equation model developed by Krishnamurti. Two 36-hour forecasts 
are made. A comparison of the forecasts is made at three-hour 
intervals up to 24 hours and at 36 hours. Both forecasts are compared 




to the reported winds at 12-hour intervals. 
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1. INTRODUCTION 



It is now evident that the use of pure quasi-geostrophic equations 
of motion in numerical forecasting has reached a plateau. There is 
little evidence that they will be much more accurate as far as prognosis 
is concerned. W. Blumen (1) has given evidence of shortcomings in 
the quasi-goestrophic approximations. These are due to excessive 
non-linear "adjustment times" before pure geostrophic balance is 
reached. This balance is a requirement in all quasi-geostrophic models 
but has special significance in the non-linear baroclinic models. 

Charney (2) implies that although the quasi-goestrophic models, by 
eliminating the high frequency noise disturbances, have given us valuable 
insight to the physical properties of meteorological significant waves, 
they are only the beginning of a new era in numerical forecasting. He 
mentions the inherent errors of the geostrophic approximation and that 
further improvements in the models would only lead to highly complex 
methods which would tend to destroy the basic idea of the geostrophic 
approximation -- its simplicity. Charney suggests that the primitive 
Eularian equations as first proposed by Richardson are the next step 
in successful numerical forecasting. 

It is known that extremely accurate initial wind conditions are 
needed in order to successfully determine accelerations and divergence 
from the primitive equations; acceleration and divergence are an order 
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of magnitude smaller than the measurable quantities in the atmosphere 
(3). Use of the primitive equations will then depend on our ability to 
adequately describe an initial wind field that will: 1) be a close 

representation of the actual wind, 2) be compatible with numerical 
calculation as far as stability and time incrementation is concerned, 
and 3) be obtainable from present measured data. Charney (2) develops 
the "Balance" equation which is used to describe an initial wind field 
that meets the above criteria. (See Section 2 for this development). 

Charney also considers the geostrophic wind components, u and 

O 

V , as part of the initial wind conditions in the primitive equations. 

O 

Through a numerical example he shows that large inertio-gravitational 
oscillations occur. These results show the existence and accumulation 
of large amounts of divergence, a characteristic of gravity waves but 
not of the meteorological significant waves. 

In this paper the "Balance" equation is developed. Using various 
modifications of the equation, four different non-divergent stream 
functions are determined and compared with their isotach patterns. 
Using Krishnamurti' s primitive equation model, two 36-hour forecasts 
are made, one with the initial winds being geostrophic (a modification 
of the "Balance" equation), and the other with initial winds determined 
from the unmodified "Balance" equation. The two forecasts will be 
compared at three-hour intervals up to 24 hours and at 36 hours. Both 
forecasts will be compared to reported winds at 12, 24 and 36 hours. 
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2. DEVELOPMENT OF THE NON-DIVERGENT WIND EQUATIONS 



There were two approaches to the development of the non-divergent 
wind equation, henceforth called the "Balance” equation. Charney*s 
approach (2) was that of necessity. He found that something better than 
the geostrophic approximation was needed in order to determine the 
initial wind fields for use in the primitive equations. Examining the 
actual divergence in the atmosphere, he noted that horizontal divergences 
were very small. Therefore, he said, the horizontal wind components 
can be described accurately as partial derivitives of a non-divergent 
stream function, \jj , as follows: 



By placing these values into the divergence equation, then eliminating 
insignificant terms, Charney arrives at the "Balance" equation. He 
implies that since this equation is a representation of non-divergent 
winds, and since the inertiogravitational winds are characterized by 
relatively large divergences, then the equation, if satisfied, will 
prohibit the occurrance of gravity waves. He mentions collaborators 
who have shown the same results and proved his implications. 

A more sophisticated approach to the development is that of 
Thompson (5). His intention was to re-examine the filtering problem 
of the geostrophic approximation and to see if it could be improved by 
providing the filtering (elimination of gravity waves) with less 
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restrictions. To do this it was necessary to determine the feature of 



the equation of motion that generates the gravity waves, and then to 
eliminate this feature with some kind of filtering approximation. 
Thompson accomplishes this by obtaining the frequency equation of the 
linearized (by perturbation method) divergence equation. He determines 
that the roots of the frequency equation which correspond to gravity 
waves, are those due to the total derivative of the divergence in the 
divergence equation. By removing this total derivative there is 
eliminated the possibility of having any gravity waves when the modified 
divergence equation is satisfied. This modified equation is the ’'Balance” 
equation as determined by Charney. 

From Helmholtz, velocity can be expressed as the sum of non- 
diver gent and divergent parts: 

V = kxV^ + vx . (2-1) 

Starting with the horizontal equation of motion 



dt 



/H 



- f kxv , 



( 2 - 2 ) 



the horizontal divergence equation is obtained: 



<3 V- V ^ u^(V- V) + 
"5~t 3 X 



v ^(V» V) 

h V 



+ ^ u • ^ ^ » 

a X 



dv _ 
J~y~ 




(2-3) 



Removing the total derivative of horizontal divergence from (2-3) leaves, 

= . (2-4) 

d ^ o y ^ 

Substituting the following expressions into equation (2-4): 



V = kxy ^ 




(2-5) 

O X 
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the ’’Balance" equation of the non-divergent wind is obtained: 









(2-6) 



L“T^ 'TTJ ‘ 

This equation defines the non-divergent wind from a height field, z, 
and eliminates inertio^ravitational waves as long as it is satisfied. 

Using equation (2-6) Krishnamurti (4) develops expressions for 
the following four stream functions: 

(1) the geostrophic stream function by dropping the non-linear 

term on the right side of (2-6) and treating the coriolis parameter, f, 



as a constant, giving 

r2.,, _ 1 „2 






or 

(2) f j^, the linear balanced stream function by dropping the non- 
linear term and treating f as a variable, giving 

- T +tV^$ ■ 

(3) the semi-geostrophic balanced stream function by 
evaluating the non-linear term using geostrophic wind components 



giving 



VY =-T <Vf*Vf) + TV"^ 



(4) complete balanced stream function by making no 

changes or assumptions in the ’’Balance’’ equation. 

The numerical considerations of convergence in solving for the com- 
plete balanced stream function are discussed by Krishnamurti. (2 -6). 
Generally, the requirement for a convergent solution is that the equation 
be elliptic. This requirement is met if: 

vY>l • 
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3, PROCEDURES 



Using 200-mb data for 1200Z, 13 April 1964, an accurate 
geopotential height field (50 meter interval) was plotted within the 
following area: 25N to 60N, TOW to 135W. These contours represent 
the geostrophic stream function, y/g* Krishnamurti (4) uses this 
height field to compute the stream functions, 

Figure 1 is the geopotential height field. Figures 2, 3 and 4 show the 
stream functions, y/j^, y/jNjj y^N2' respectively, plotted with the 

height contours on a polar stereographic projection. Stream function 
interval is 50 X 10^ m^ sec"^. 

Taking the expression for the non-divergent wind: 

V = k xvy/ , 

wind speeds were computed at over 200 points in the fields of 

and ^N2. and isotachs were drawn, Isotachs of the geostrophic 
wind were drawn to speeds computed from: 

V = k X ^ . 

Figures 5 through 8 are the isotach fields of the four stream functions. 
Three graphs of wind speed versus latitude were drawn. The 
first is along the line through the isotach maximums, labled line "A"; 
the second along the 95th meridian, labeled line the third along 

the 115th meridian, labeled line *'C", These graphs are Figures 10, 

11 and 12 respectively. 
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Reported winds are shown in Figure 9. Figure 13 is a table 
comparing the reported winds to each of the four computed winds, 

4 . DISCUSSION OF WIND FIELD COMPARISONS 

Examining the streamline patterns, Figures 2, 3 and 4, we see 
that definite cross-contour flow exists. Generally the flow is to lower 
heights upstream of close contour spacing and to higher heights where 
contours spread. The linear balanced stream function, being closest 
to the geostrophic wind in computational simplicity, exhibits the least 
cross-contour flow. The complete balanced stream function shows 
the most cross-contour flow and therefore corresponds more closely 
to statistical observations of this phenomenon (6). The flow exhibited 
by the non-linear stream functions is a direct manifestation of the 
theories on geostrophic adjustment (7). 

Only the complete balanced stream function shows strong anti- 
cyclonic flow over New Orleans, This is the area of the sub-tropical 
jet stream. The anti-cyclonic flow seems to be excessive but may be 
due to boundary values in computations or a lack of data. 

It must be remembered that the streamlines represent the non- 
divergent wind. The divergent parts of the wind are small. Krishnamurti 
says that calculations for this map time give divergent winds that are 
at no place greater than five percent of the total wind. 

Streamline spacing corresponds well to contour spacing (geostrophic 
streamlines). This is reflected in the isotach patterns of Figures 5 
through 8, and the three graphs, Figures 10, 11 and 12, The linear 
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balanced winds are close to geostrophic winds. Krishnamurti has 
shown (4) that at low levels where strong cyclonic and closed circula- 
tions exist, the linear balanced winds are far too large. At 200-mb 
the geostrophic wind maximum in the trough will approach the complete 
balanced wind when corrected for curvature. The semi-geostrophic 
balanced winds are close to the complete balanced winds. However, 
only the complete balanced winds give a good representation of the sub- 
tropcial jet. 

Figure 13 shows a numerical comparison of the computed and 
reported winds. The overall mean algebraic difference indicates that 
the balanced winds are closer to the reported winds than the geostrophic 
winds are. A comparison of this type at the 200-mb level must be 
viewed with caution. It is known that high values of vertical wind shear 
exist at this level. Small errors in balloon height can give unrepresenta 
tive reported winds* Looking at the isotach patterns of the computed 
winds, we also see large horizontal shears. In areas of high winds and 
large shears, the downwind position of the balloon from the reporting 
station could be in a sector of the computed wind isotach field that 
contains a velocity range of up to 20 knots. It would be very difficult 
to include the turning of the wind with height to get an accurate balloon 
position in order to compare the reported and computed wind at that 
point. It might be advantageous to make a smooth isotach analysis of 
the reported winds and use this to compare the computed winds. There 
still will exist, however, the problem of balloon location and which 
winds have errors. 
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Lacey (9) compared the geostrophic winds and the non-divergent 
complete balanced winds with the reported winds at 850 -mb, 700-mb, 
500-mb and 300-mb over the United States, His comparison was made 
for both direction and speed. He found that on the average the com- 
plete balanced winds were much closer to the reported winds in both 
direction and speed when the complete balanced wind and the geostrophic 
wind differed by more than five degrees or five knots. He also 
determined that adding the divergent wdnds to the non-divergent complete 
balanced winds produced little change in the non-divergent wind pattern. 
The divergent winds were about two to three percent of the total wind. 

From the above discussion it can be concluded that the complete 
balanced wind is the most accurate description of the true wind. The 
question, “are the complete balanced winds accurate enough for use 
as initial wind conditions in the primitive equations? “ will be discussed 
in the next section. 

5. THE NUMERICAL FORECASTS 

In the Introduction it was mentioned that Charney proposed 
balanced winds as sufficient initial conditions for numerical prognosis 
using the primitive equations. His main point was that the geostrophic 
wind did not filter the “noise" while the balanced winds did. He used a 
simple model to prove this point. However, Phillips (8) has shown that 
because Charney^ s model was too simple, he overlooked the fact that 
using only the non-divergent winds as initial conditions (V* V = 0 at 
t = 0), the forecast using an actual baroclinic atmospheric model 
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would still contain considerable ’’noise”. Phillips claims that the noise 
generated by the restriction to non-divergence could be eliminated by 
including the divergence obtained from geostrophic approximations in 
the initial data. 

Krishnamurti has developed a primitive equation model in which 
he has incorporated the ideas of Charney, Phillips, Thompson and 
others. His initial wind field is a combination of the non-divergent 
wind determined from the ’’Balance” equation and the divergent wind 
obtained by computing divergence. Any combination of computing 
methods to calculate balanced winds and vertical motions (thus, 
divergence) can be used to determine the initial conditions. He has 
done extensive work on partitioning the atmosphere by using his model 
with initial conditions determined as follows: 

a. The non-divergent stream function, using the 

complete non-linear balance equation. 

b. Vertical velocity expressed in a complicated omega equation 

involving divergent velocity potential, X ^ 

c. The divergent velocity potential, X * expressed in the mani- 
festation of the continuity equation: 




From the initial conditions obtained in the above manner, a 36-hour 
forecast was obtained. 

A much simpler procedure would be to use the quasi-geostrophic 

approximation in all phases of calculating the initial wind conditions. 

Though this approximation itself has proven to be inadequate if used 
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solely as the initial wind condition, one might ask that if it were com- 
bined with quasi-geostrophic divergence as suggested by Phillips, 
would it still be inadequate? A second forecast was made using 
Krishnamurti* s model with the initial conditions determined as follows: 

a. The non-divergent stream function, \j/ , determined from the 
geostrophic approximation, 

b. Vertical velocity expressed in the quasi-geostrophic omega 
equation. 



c. The continuity equation expressed as: 



g 



d P 



In both forecasts, grid distance was: 

y = 2-1/2 degrees latitude, 

X = y sine latitude. 

The area covered was from 55W to 135W and from 25N to 60N. Time 
steps were ten minutes. North and south boundary conditions were 
fixed. The eastern boundary was extended eastward six grid distances 
and combined cyclically with a w^estern boundary extended six grid 
distances westward. Data used were the height field from OOOOZ, 

12 April 1964. Figures 14, 15 and 16 show the contours, geostrophic 
wind isotachs, and reported wind isotachs, respectively. 



6. DISCUSSION OF RESULTS OF THE FORECASTS 



Figures 17 through 25 show the complete balanced and quasi-geo- 
strophic component u and v fields at three-hour intervals up to 24 hours and 
at 36 hours, respectively. These figures yield the following observations: 
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a. At three hours (Figure 17) there are already considerable 
differences between the balanced and geostrophic components. 

In the u field at 55N, 75 W, there exists a negative 5 m-sec"^ 
center in the geostrophic field, but a positive 15 m-sec ^ in the balanced 
field. The elongated center on both over the northwestern United States 
if 15 m-sec"^ larger in the balanced field. The 50 m-sec" ^ center off 
Washington in the balanced field is only 20 m-sec" ^ in the geostrophic 
field and is displaced 15 to 20 degrees west. The sub-tropical jet over 
Florida is expressed well in the balanced field but is displaced south 
in the geostrophic field. Many irregularities appear in the geostrophic 
field while the balanced field exhibits a smooth, realistic pattern. 

In the V fields there is much more agreement. This should be 
expected because of the fixed boundary values. We note that the 
positive area centered at 105W on both fields is much smaller in area 
and strength in the balanced field. The opposite is true for the negative 
area at 125W. A look at Figures 15 and 16 reveal that the balanced 
field is more realistic. 

b. At six hours (Figure 18) the differences noted at three hours 
between balanced and geostrophic components are even more pronounced. 
There are now strong indications of inertio- gravitational oscillations 

in the geostrophic field. Lookingat the u fields we notice that a large 
high latitude negative center extends the length of the chart in the geo- 
strophic field, while the balanced field shows only a decrease in the 
positive zonal flow in the same area. The similarity over northwestern 



20 



United States which existed at three hours has almost disappeared even 
though the balanced field shows little change during the three hours. 

The balanced field of the v component shows little change in three hours, 
but the geostrophic v field shows a great increase of negative components 
at higher latitudes. The overall similarity that existed at three hours 
is disappearing at six hours. 

c. At nine hours (Figure 19) the balanced components are 
behaving well. There is a slight easterly progression of the u maximum 
at 130W, 45N, and the ridge line in the v field along the west coast has 
easterly movement. In the u field we see the sub-tropical jet maximum 
moving eastward. 

The geostrophic u field has become strongly negative above 
SON. The indication of a sub-tropical jet is becoming distorted, and 
it is becoming difficult to distinguish the developing trough. The v 

field now shows good continuity with the past three hours and is becom- 

1 

ing similar to the balanced v field. 

d. The balanced components show good continuity at 12 hours 
(Figure 20). The geostrophic u field is now rapidly losing the high 
latitude negative center that existed at nine hours. The geostrophic v 
field retains continuity but there is evidence of oscillations. The 
negative center at 125W is lagging the balanced negative center by 
about ten degrees. 

e. The balanced components at 15 hours (Figure 21) still have 
good continuity and retain pattern shape. The geostrophic v field seems 



21 



to have settled down with a pattern development similar to the balanced 
field. The geostrophic u field has lost its negative values but still 
shows oscillations. This field has only half the positive velocity as 
the balanced u field, and no marked centers exist. 

f. For 18 hours (Figure 22) and 21 hours (Figure 23) the same 
trends exist. The geostrophic u field has again established a negative 
area above SON while the rest of the field bears little resemblance to 
the balanced one. 

g. At 24 hours (Figure 24) the geostrophic u field has the negative 
area in high latitudes with no indication of a jet stream anywhere. 

The pattern slightly resembles the balanced pattern but the values are 
all much lower. The balanced u field appears to be developing strong 
positive values in the lower left corner. This may be due to the cyclic 
feed through from the lower right maximum which existed at the start 
of the forecast. The v fields show pattern similarity but the geostrophic 
field has a great predominance of positive values; the negative areas 
having small values. 

h. The 36-hour prognostic u and v fields (Figure 25) show a 
pattern breakdown. The balanced u field has strong distortion from 
cyclic feed through, but maintains a resemblance of the synoptic 
situation in its pattern shape. The geostrophic u field shows strong 
distortion with a negative area undercutting a positive one. The balanced 
V field has distortion in the western half with weak centers indicating a 
very weak ridge. The geostrophic v field shows a counter-clockwise 
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rotation of the prevailing pattern in the western half, with a very 
strong negative center just east of Florida. 

It can be concluded from the above discussion that the quasi- 
geostrophic approximation, even by including divergence, is unstable 
and gives rise to large inertiogravitational oscillations when used for 
initial conditions in the primitive equations. From the limitations and 
boundary values applied in this model, the conclusion is most evident 
in the u component field. 

A further comparison of the prognoses was made by comparing 
the reported wind speeds to the prognostic winds. Prognostic winds 
were obtained by vector addition of the u and v fields. Figures 26, 27 
and 28 show the reported winds and the prognostic balanced and geo- 
strophic winds at 12, 24 and 36 hours, respectively. The table of 
Figure 29 shows definitely that the balanced winds are much better than 
the geostrophic winds for the 12 and 24-hour prognoses. At 36 hours 
too much distortion and cyclic feed through exist to make a reliable 
comparison, for the geostrophic prognosis appears to be better. 

The major conclusion in this numercial forecast comparison is 
that the balanced winds give the prognosis a wind pattern that resembles 
the synoptic situation. The geostrophic winds gave rise to considerable 
large-scale ’’noise” and unrealistic patterns. 
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